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Abstract 

Imaginary-time Nambu-Bethe-Salpeter (NBS) wave function is introduced to 
extend our previous approach for hadron-hadron interactions on the lattice. 
Scattering states of hadrons with different energies encoded in the NBS wave- 
function are utilized to extract non-local hadron-hadron potential. "The 
ground state saturation", which is commonly used in lattice QCD but is 
hard to be achieved for multi-baryons, is not required. We demonstrate that 
the present method works efficiently for the nucleon-nucleon interaction (the 
potential and the phase shift) in the 1 5'o channel. 
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1. Introduction 

Euclidean correlation functions are dominated by contributions from the 
corresponding lowest-energy states at sufficiently large time separation t. 
This property, called the ground state saturation, is heavily used in lattice 
QCD to extract various hadronic observables such as masses, decay constants 
and other matrix elements. The ground state saturation, however, is difficult 
to be achieved for multi-baryon systems. For example, the signal-to-noise 
ratio for the correlation of n-nucleons reads ll|: 



] ^ e -n(m N ~3m 7r /2)t /-j\ 

/ n 

where m^ and m n are the nucleon mass and the pion mass, respectively. 
This relative enhancement of statistical noise at large t for m^ — 3m 7r /2 > 
is a common problem for baryonic systems, even for a single baryon (n — 1). 
In addition, there exists another problem for the multi-hadrons at large t 
caused by the small splitting between the ground and the 1st excited states 
for large volume. For example, in the nucleon-nucleon (NN) system, the 
minimum splitting is estimated as 

A£^ = J-^, (2) 

where L is a spatial extension of the lattice. If L ~ 6 fm and m^ ~ 1 GeV, 
we have AE ~ 43 MeV ~ 1/(4.6 fm). The ground state saturation requires 
t ^> (AE)^ 1 ~ 4.6 fm, which corresponds to t/a 3> 46 for the lattice spacing 
a ~ 0.1 fm. It is very difficult to extract signals at such large t due to the bad 
behavior of statistical noise in Eq. ([T|). To avoid these problems, techniques 
such as the use of improved operators and/or the diagonalization of matrix 
correlation functions [2| have been employed. U 

Recently, a novel method to derive hadron-hadron interactions from lat- 



tice QCD was developed by HAL QCD collaboration [5|-[TJJ, where the 
Nambu-Bethe-Salpeter (NBS) wave function is utilized to extract the hadron- 
hadron potentials. Since its extraction from hadronic correlation functions 



^or recent applications of these methods to multi-baryons, see e.g. Ref. [3, | 



relies on the ground state saturation, the problems mentioned above may ex- 
ist in principle. In this paper, we introduce time dependent Schrodinger-like 
equation which can utilize the information of moderate t and can avoid the 
problem of ground state saturation. This allows one to derive the hadron- 
hadron potentials as defined in the original HAL QCD method with less 
systematic errors. The key observation is that the scattering states with 
"different" energies on the lattice are governed by the "same" non-local po- 
tential U(r, r 1 ). The present time-dependent method has already been ap- 
plied successfully to the baryon-baryon (BB) potentials in the flavor-SU(3) 



limit [12|, |13j. In the following, we give a full theoretical account of this 
imaginary-time HAL QCD method by taking the NN scattering with (2+1)- 
flavor lattice QCD as a concrete example. 

2. Potential with Ground State Saturation 

In the original HAL QCD method, the "time-independent" NBS wave 
function was shown to satisfy the following "time-independent" Schrodinger 
equation [7j]: 



(k 2 /m N 



Ho) M^ = / d 3 r'U{f, f)^(0, ( 3 ) 



where H = —A/m^ with m^ being the nucleon mass. The potential U(r, f 



is non-local but independent on k [7J, ll4|. The equal-time NBS wave function 
is given by 

iP % (x-y) = (0\N(x)N(y)\N(k)N(-k);in), (4) 

where |0) and \N(k)N(— k);in) denote the vacuum and a two-nucleon state 
with an asymptotic momentum k, respectively, and N(x) denotes a compos- 
ite interpolating field for the nucleon. For N(x) being local, the reduction 
formula for local composite operators [15| can be used to establish the re- 
lation between the NBS wave function and the S-matrix: The asymptotic 
behavior of the NBS wave function at long distance reads [7|, [l6|, [17 



, ,.. iS(k) sm( kr + 5(k)) ,. 

tl)%{r) = e im — — + ■ • • , (5) 

where 5(k) denotes the (scattering) phase of the S-matrix. Therefore the 
NN potential in Eq.Q gives correct phase shift 5(k) for all k in the elastic 
region E < E t ^ = 2m n + m n , by construction. 



In lattice QCD calculations, NBS wave functions are extracted from the 
NN correlation function using the ground state saturation as 

C NN (x-y;t) = ^Y,(°\ T {N(x + ?,t)N(y + r,t)-j(0)}\0} 

r 

= J2^n(x - y) ■ a n e' Ent 

-* xjj (x - y)a e- Eot (t->oo), (6) 

where J7(0) is a two-nucleon source located at t — 0, V denotes the spa- 
tial volume, ip n (x — y) = (0\N(x)N(y)\n) denotes an NBS wave function 
for an intermediate state \n) with the energy E n , and a n = (n|j7"(0)|0). 
In the spin-singlet sector, for example, the central potential Vc(r) in the 
leading order of the velocity expansion, U(f,f f ) = V(f, V r )5 3 (r — f") = 
{Vc(r) + 0(V 2 )} <5 3 (f — r"), is given by 

ir f \ k H C NN (r,i) 

Vcir) = hm — - — 7^-^, (7) 

where k denotes the "asymptotic momentum" for the ground-state. The 
ground state saturation is crucial here to extract the potential. 

3. Potential without Ground State Saturation 

In this section, we propose an alternative derivation of potential with- 
out using the ground state saturation. For this purpose, we consider the 
normalized NN correlation function 

R(t,f)=C NN (t,f)/(e- mNt ) 2 . (8) 

We here assume that t is moderately large such that elastic contributions (at 
E < E th = Ira N + m n ) dominate Cjvjv(i, 0- 
As before, we write 

R(t,r) = 5>je(0 • «^xp (-tAW(k)) , (9) 



where AW(k) = 2^/m 2 N + k 2 - 2m N and a^ = (N (k) N(-k);in\ J (0)\0). 
From an identity AW(k) = — A , it is easy to see 

= e {*» + f - ^} wo • «i«p (-««; 

where C/ is the integration kernel associated with the non-local potential 
U(f, f). Thanks to Eq.flB]), k 2 /m^ in the first line can be replaced by H + U 
in the second line. We then arrive at the "time-dependent" Schrodinger-like 
equation, 

{i^ - m ~ H °) ^0 = J * 3 r>u(r,n m ,n, {u) 

which shows that the same potential U(r, r 1 ) as defined in Eq. ([3]) can be 
obtained from R(t, r). An advantage of this method is that the ground state 
saturation (or more generally a single state saturation) is not required for 
R{t,r) to satisfy Eq.(HI]). 

Several comments are in order here: 
(i) For the present method to work, t has to be large enough such that elastic 
contributions at E < E t ^ = 1m n + m n dominate R(t,r). Note that such t 
is much smaller than that required for the ground state saturation. This is 
especially so for the large volume, since typical size of gaps between energy 
eigenvalues shrinks as 0(1/L 2 ). While it becomes more and more difficult 
to achieve the ground state saturation for larger volume, the requirement of 
the elastic dominance E < E t t is less sensitive to the volume size, 
(ii) In the leading order of the velocity expansion, Eq.( 1TT|) leads to a gener- 
alization of Eq. © 

H R(t,r) (d/dt)R(t,f) 1 (d/dt) 2 R(t,r) 

CU " R(t,r) R(t,f) Am N R(t,r) [ ' 

for the spin-singlet sector. We can also include the higher order terms of the 
velocity expansion as discussed in Ref . (3, |l0] ■ Convergence of the velocity 
expansion in the original HAL QCD method can be examined by comparing 



local potentials at two different energies as discussed in Ref. [10[. Equiva- 
lently, in the present method, the convergence can be examined by comparing 
the local potentials obtained for different t's. 

(iii) D = j^^ - | in Eq.(HI]) plays a role of k 2 /m N in Eq.©. For 
AE ■ t ^> 1 where the ground state saturation is achieved, Eq.flTT]) reduces to 
Eq.(in]). We therefore can regard the "time-dependent" Schrodinger-like equa- 
tion as an extension of the time- independent Schrodinger equation (Eq.©). 



m N 



(iv) In the non-relativistic limit where AW(k) = 2ym 2 N + k 2 — 2m n 
"time-dependent" Schrodinger-like equation leads to 

|-^ - H \ R(t,f) = J ' d 3 r'U(ry)R(ty). (13) 

Therefore, the 2nd derivative term of t in Eq. ( TTTT) corresponds to the rela- 
tivists effect. 



4. Numerical Results 

To test the present method, we employ (2+l)-flavor QCD gauge configu- 
rations generated by PACS-CS collaboration [18| on 32 3 x 64 lattice with the 
RG improved Iwasaki gauge action at f3 = 1.9 and the non-perturbatively 
0(a) improved Wilson quark action at (K ud) n s ) = (0.13700,0.13640) and 
Cgw = 1-715. This parameter set corresponds to the lattice spacing a ~ 0.091 
fm (a- 1 = 2.176(31) GeV), the spatial extent L = 32a ~ 2.90 fm, m n ~ 701 
MeV and m N ~ 1583 MeV. 

The periodic boundary condition is used for spatial directions, while the 
Dirichlet boundary condition is taken for the temporal direction at £dbc = 
32a and —32a, to avoid opposite propagations of two nucleons in tempo- 
ral direction, i.e, one propagates forward and the other propagates back- 
ward. From time-reversal and charge conjugation symmetries, we can aver- 
age over forward propagation at t > and backward propagation at t < 
with a source at t — 0. By temporally shifting gauge configurations, 21 
source points are used per one configuration and 390 gauge configurations 
are employed in total. Statistical errors are estimated by the Jackknife 
method with a bin size of 10 configurations. Composite operators for the 
proton and the neutron are taken to be p(x) = e a b c (u^^C^d^x)) u c (x) 
and n(x) = e a b c {u^(x)C'-)^di ) (x)) d c (x), respectively. In our actual calcula- 
tion, we replace e _mjv< in Eq.([8]) by the single- nucleon correlator C^{t) = 
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^2$(0\T[N(x)N(0)]\0) to suppress statistical noises of R(t,r). This replace- 
ment is permitted as long as the ground state saturation for Cat(£) is achieved. 
Time derivatives in Eq.f JT2|) are evaluated after applying the polynomial in- 
terpolation of degree 5 to R(t,f). 

In order to see how excited states of the two nucleons affect the final NN 
potential, we introduce a source function with a real parameter a as 

f(x, y,z) = l + a (cos(2ttx/L) + cos(2iry/L) + cos(2nz/L)) , (14) 

which reduces to the wall source at a = 0. Two-nucleon source is then 
defined by J(f) = p(f)_- n(f),_ where p(f) = e abc (u a (f)C'y 5 d b (f)) u c (f) and 
n(f)_= e abc (u a (f)Cr/ 5 d b (f)) d c (f) with u(f) = J2x U (^)f(^) and d U) = 

Figure [U (left) shows C NN (r,t) at t = 9 for a = 0.00, 0.08 and 0.16. If 
the ground state saturation were achieved, results with different values of a 
should be the same up to overall normalizations. Figured] (left) reveals that 
contamination from the excited states is non-negligible at t — 9. As shown 
in Figure [T^right), the contamination is transferred to the a dependence of 
[H C NN (f,t)}/C NN (f,t). 

Figure [2](left) shows Vc(r) obtained from our present method Eq. ffT2|) for 
three values of a. The a dependence seen in Figured^ right) disappears within 
statistical errors. Three contributions to Vc(r) in Eq. (ll2p are separately 
shown in Figure fright) for a = 0. We observe that the first term of Eq. (TT2|) 
(the red points) determines the main trend, while the second term (the blue 
points) gives an important correction. The third term (the green points), on 
the other hand, is negligible, showing that the non-relativistic approximation 
AW(k) ~ k 2 /rriN works well in this case. Note that the ^-dependence of the 

second term in Eq.( TT2|) . — ( ' ^ i, ' r > = ° g a/ ? is a useful measure of the 

departure from the ground state saturation. 

5. Scattering Length and Scattering Phase Shift 

We now calculate the A^A^ scattering phase shift, by solving the Schrodinger 
equation with the potential Vc(r) in the infinite volume. For this purpose, 
as shown in Figure EJleft), the central potential Vc(r) is fitted with multi- 
Gaussian functions as g{r) = J2n=i BB Vn ' ex P(~ ^n^ 2 ), where V n and u n (> 0) 
are used as fit parameters, A gauss denotes the number of Gaussian functions. 



We then solve the Schrodinger equation in x Sq channel 19 



Figure fright) shows the scattering phase 5(k) extracted from the long 
distance behavior of the solution ipk{ r ), together with the experimental data 
for comparison. Qualitative feature of the phase shift as a function of k is well 
reproduced, though the strength is weaker, most likely due to the heavy pion 
mass (m-ir ~ 701 MeV) in this calculation. In fact, the recent 3-flavor QCD 
simulations show that the NN phase shift approaches toward the physical 



value as the quark mass decreases [13|. The scattering length for m w ~ 701 



MeV in the present method, calculated from the derivative of the scattering 
phase shift at E\ & y, = 0, leads to a( 1 S'o) = lim^o tan £(&)/& = 1.6 ±1.1 fm. 

6. Non-local potential 

By extending the present method further, one may directly extract the 
non-local potential. Let us introduce a more general A^ A" correlation function, 

1 ^ (0\T[N(x + r t)N(x, t)N(x' + f, 0)N(x>, 0)]|Q> 

Rr,r>(t) = —^ {e -m N ty ' ( 15 ) 

x,x' 

which is shown to satisfy 

m s {^S-li- H °} m = u - m - (16) 

Here the matrix indices r, r" for K(t), R(t), U and a necessary integration over 
spatial coordinates are implicit. The non-local potential is then extracted as 
U = K(t) ■ i? _1 (t), where an approximated inverse of the hermitian operator 
R(t) is defined by -R~ 1 (t) = J2\ n <t)=£o ^n(^)~ 1 | r M)( r M|- Here \ n (t) and \n,t) 
are an eigenvalue of R(t) and its eigenvector, respectively. Zero eigenvalues 
are removed in the summation. Note that U which satisfies Eq. flTB"]) is not 
unique, since U = U + J2\„(t)=o c n\n,t)(n,t\ also satisfies the same equation 
for arbitrary {c n }. This is related to the fact that the zero-modes or nearly 
zero-modes are associated with states above the inelastic threshold. 

7. Summary and Concluding Remarks 

A method to extract hadron-hadron interactions by generalizing the origi- 
nal HAL QCD method is proposed. We derived "time-dependent" Schrodinger- 
like equation, a second order differential equation in t, which enables us to 
construct A^A^ potentials without assuming the ground state saturation in 

8 



the elastic region E < E t ^ = 2m^ + m^. We have shown that this method 
works well for extracting the central NN potential in the 1 Sq channel: Iden- 
tical potential is obtained for different source-operators within the statistical 
error. Also, resultant NN potential, the scattering phase shift and the scat- 
tering length are much improved both qualitatively and quantitatively from 
those obtained by assuming ground state saturation. 

While we have considered the system only in the elastic region so far, 



an extension to the inelastic region is also possible. In Ref. [21[, it has been 
shown that one can define and extract the hadronic potentials above inelastic 
threshold. Together with the framework proposed in this paper, this could 
provide a novel prescription to solve the S-matrix in QCD. 
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Figure 1: (left) C NN (r,t) at t = 9. (right) -H C NN {r,t)/C NN {r,t) at £ = 9. 
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Figure 2: (left) Central potentials obtained by our new method Eq. (fT2|) at £ = 9 for three 
values of a. (right) Three contributions to Vc(r) in Eq. ([T2"l) at t = 9 for a = 0. 
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Figure 3: (left) The multi-Gaussian fit of the central potential Vc(r) with iVcauss = 5 
for a = at t = 9. (right) The scattering phase in So channel in the laboratory frame 
obtained from the lattice NN potential, together with experimental data[20(. 
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